Mechanistic simulations of population dynamics

We simulated antagonistic population dynamics using a discrete Lotka-Volterra competition model, where each population’s growth rate is influenced by the fluctuations in a competitor population. Pairs of populations were simulated according to:

\[\begin{equation} N_i(t+1) = N_i(t) + \lambda_i * N_i(t) \end{equation}\]

where \(N_i\) is the size of population \(i\), \(t\) is each time step, and \(\lambda\) is the true growth rate of population \(i\). The true growth rate (\(\lambda\)) is determined by:

\[\begin{equation} \lambda_i = r*(1-\frac{(N_i(t) + \alpha_{ji}N_j(t))}{K_i}) \end{equation}\]

where \(r\) is population \(i\)’s maximum growth rate, \(\alpha_{ji}\) is the interaction effect of population \(j\) on population \(i\), \(N_j\) is the size of population \(j\), and \(K_i\) is the carrying capacity of population \(i\).

To test index sensitivity to error, we introduced process and observation error into the generated population time series.

Table 1: Constant simulation parameters
Value
Number of population pairs 10.0
Time steps 10.0
Initial size (N0) 100.0
Maximum growth rate (r) 1.5

Biodiversity change scenarios

Direction of change

To test index sensitivity to direction of change, we generated three scenarios of population-level biodiversity change.

We varied carrying capacity through time to drive populations to:

    1. decline in response to decreasing carrying capacity, representing situations where habitats are lost or altered,
    1. remain stable in response to constant carrying capacity, where habitats are left untouched, or
    1. increase in response to increasing carrying capacity, where habitat quality and/or availability increases through time.

In each scenario, carrying capacity begins at \(K = 100\) and each population’s initial size is exactly at this carrying capacity (\(N_0 = 100\)).

Carrying capacity trends used to generate the biodiversity change scenarios.

Figure 1: Carrying capacity trends used to generate the biodiversity change scenarios.

Covariance

To test whether the LPI is sensitive to lags in the covariance between populations, we then generated two levels of lag in the interaction effect dictating how populations covary through time:

  • The first level was lag-1, where the interaction effect of populations \(i\) on populations \(j\) was dependent on the abundance of populations \(i\) at the previous time step (\(t-1\)).
  • The second level was lag-2, where this interaction effect depended instead on the abundance of populations \(i\) from two time steps ago (\(t-2\)).

Measuring index sensitivity

For each biodiversity change scenario, index sensitivity was measured as the accuracy and precision of the calculated Living Planet Index when compared to the true biodiversity change trend. This true LPI trend was obtained by calculating the index values through time based on a population generated from the same parameters we used to generate populations for each scenario, without any added process error.

We measured accuracy as the difference between the calculated LPI and true LPI trend at each time step, expressed as a percentage of the true LPI value at each step to allow for comparison between scenarios:

\[accuracy = \frac{LPI_{calculated} - LPI_{true}}{LPI_{true}} * 100% \]

As a first metric of precision, we evaluated whether the true LPI trend fell within the 95% confidence interval around the calculated LPI trend obtained from bootstrap resampling at each time step. As a second precision metric, we measured precision as the difference between the expected uncertainty interval and the width of the 95% confidence interval expressed as a percentage of the calculated LPI value at each step to allow comparison between scenarios. The expected uncertainty interval is determined based on the process and observation errors introduced into the simulated population trends. I calculated each population’s growth rate with and without introduced error, and took the difference between the two to get a measure of the total introduced uncertainty in the growth rates. Then, I converted the difference between these growth rates to the same scale as the LPI:

\[ precision = \frac{LPI_{CIwidth} - LPI_{uncertainty}}{LPI_{calculated}} * 100% \]


Is the LPI sensitive to direction of change?

Table 2: Simulation parameters
i j
Interaction effect (alpha) 0 0
Observation error +/- 5 +/- 5
Process error +/- 0, 0.1, 0.2 +/- 0, 0.1, 0.2
Lag 0 0

Process error = 0

True and calculated LPI trends for each direction of biodiversity change.

Figure 2: True and calculated LPI trends for each direction of biodiversity change.

Process error = 0.1

True and calculated LPI trends for each direction of biodiversity change.

Figure 3: True and calculated LPI trends for each direction of biodiversity change.

Process error = 0.2

True and calculated LPI trends for each direction of biodiversity change.

Figure 4: True and calculated LPI trends for each direction of biodiversity change.

Summary

Process error = 0

LPI sensitivity to direction of change.

Figure 5: LPI sensitivity to direction of change.

Process error = 0.1

LPI sensitivity to direction of change.

Figure 6: LPI sensitivity to direction of change.

Process error = 0.2

LPI sensitivity to direction of change.

Figure 7: LPI sensitivity to direction of change.

The LPI is sensitive to direction of change:
* declining trends are overestimated by the LPI 
* declining trends are estimated with lower precision

Is the LPI sensitive to covariance?

To test whether the LPI is sensitive to covariance between populations, we generated biodiversity change scenarios with two types of interaction effects (\(\alpha\)) between populations:

  • For the first scenario, we set a negative interaction effect to represent an antagonistic interaction between populations \(i\) and \(j\), such as a competitive or predator-prey relationship, which causes negative covariance.
  • The second scenario was based on a positive interaction effect of the same magnitude(s) to replicate a positive interaction between populations \(i\) and \(j\), and therefore positive covariance.

Each type of interaction was tested at two interaction strengths. All interaction effects were immediate, meaning there was no time lag in the covariance in the populations’ growth rates.

Table 3: Simulation parameters
i j
Interaction effect (alpha) 0 0
Observation error +/- 5 +/- 5
Process error +/- 0, 0.1, 0.2 +/- 0, 0.1, 0.2
Lag 0 0

Process error = 0

True and calculated LPI trends for each scenario.

Figure 8: True and calculated LPI trends for each scenario.

Process error = 0.1

True and calculated LPI trends for each scenario.

Figure 9: True and calculated LPI trends for each scenario.

Process error = 0.2

True and calculated LPI trends for each scenario.

Figure 10: True and calculated LPI trends for each scenario.

Summary

Process error = 0

LPI sensitivity to direction of change and covariance.

Figure 11: LPI sensitivity to direction of change and covariance.

Process error = 0.1

LPI sensitivity to direction of change and covariance.

Figure 12: LPI sensitivity to direction of change and covariance.

Process error = 0.2

LPI sensitivity to direction of change and covariance.

Figure 13: LPI sensitivity to direction of change and covariance.

The LPI is sensitive to covariance when trends are declining:  

* The LPI is particularly sensitive to negative covariance in declining trends, which are very overestimated (alpha = 0.1 and 0.2)  
* growing and stable trends are slightly overestimated when there is high positive covariance  
* however, this sensitivity does not change at different levels of process error   

Is the LPI sensitive to lagged covariance?

We introduced two levels of lag in the interaction effect that influences how populations covary through time:

  • The first level was lag-1, where the interaction effect of populations \(i\) on populations \(j\) was dependent on the abundance of populations \(i\) at the previous time step (\(t-1\)).
  • The second level was lag-2, where this interaction effect depended instead on the abundance of populations \(i\) from two time steps ago (\(t-2\)).

Both of these lags were introduced into positively and negatively covarying populations shown above.

Table 4: Simulation parameters
i j
Interaction effect (alpha) 0 0
Observation error +/- 5 +/- 5
Process error +/- 0, 0.1, 0.2 +/- 0, 0.1, 0.2
Lag 0, 1, 2 0, 1, 2

Positive covariance

Process error = 0

Process error = 0.1

Process error = 0.2

Negative covariance

Process error = 0

Process error = 0.1

Process error = 0.2

Summary

Process error = 0

Process error = 0.1

Process error = 0.2

The LPI is not sensitive to lags of 1 or 2 time steps.